
700 PROGRAMMING SYSTEMS 



SQUARE ROOT - MP FLOATING POINT 



QUALITY SOFTWARE 



DATE: April 1968 

ID CODE: BSH 

DRAWING: 390006 

LABEL: MSQR 

AUTHOR: JACQ 

SOURCE: SYM I 

OBJECT: Relocatable in Block "MATH" 



PURPOSE 

To extract the square root of a mid -precision floating point number and 
leave the result in the same format in the software registers MNT1, MNT2, 
MNT3. 

USAGE 

Calling Sequence 

The routine returns atL+2. 



L-l 


SMB 


MSQR 


L 


JSX 


MSQR 


L+l 


D 


ARG 


L+2 


Return 




Argument De 


scription 





The argument is a mid-precision floating point number occupying three con- 
secutive words in memory, beginning at address "ARG", as pointed out in 
the calling sequence. 

METHOD 

Let the argument X and its square root Vx be written in the floating point 
format: 



X 

N/3T 



Xm -2 



Xe 



\[M-2 J 
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Since the exponent E must be an integer: 

M = Xm ; E = —r- (for Xe even) 

M = -^2. . E = ^fti (for Xe odd) 

Except for zero, the argument is assumed to be normalized. An unnormal- 
ized argument is treated as zero. Therefore, M is within the range (1/4, 1): 

.25 ^ M ^ 1 

\Tm is extracted by the Newton-Raphson method. To limit the process to two 
iterations, a first linear approximation is made to yield 8 bits of accuracy. 
This is done by dividing the range of M into six equal intervals which are 
allocated distinct linear coefficients. 



A o 


= 


p . M + q 


= Mm + 2" V 


A l 


= 


^<VV 


= Mm + 2" 15 


A 2 


= 


l/2(f i+V 


= \|m + 2" 30 


ERROR CONDITIONS 





A negative argument is treated as zero. No error message is released. 

ACCURACY 

29 bits 

RESTRICTIONS 

Loading 

This subroutine locally references storage words and constants located in 
the "MATH POOL" module, and must therefore be loaded in the same 2k 
block as the pool. 

Other Routines 
DIVS, M. ZE, MPYS 
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Space Used 

82 words 

Timing (in cycles) 



Minimum 
Maximum 
Average 



Standard 

31 
730 
676 



Hardware Multiply &e Divide 

31 
232 
215 
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MSQR 



I 



Y.-X./2 



YES 



FETCH ARGUMENT 
X « X • 2 Xe 




NO 



1 



Y . = i X e + 1 > /2 
M 8 . C/2 



-*T»- 



_J 



Yi « (M/Y + Y ft )/2 



"1 



DEFINE INTERVAL OF 
M WITH INTEGER i 

M-l/8/i/HM 




GET LINEAR APPROX 
IMATION OF VTt : 
Y Q • P(i) . M + Q (i) 




Y m ( 1 + M)/2 



##- 



/i»Y • 2 Y « 
m 



C 



EXIT 



Y = 



Y ^ S ° 
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SIN AND COS - MP FLOATING POINT 



QUALITY SOFTWARE 

DATE: April 1968 

ID CODE: BSJ 

DRAWING: 390007 

LABEL: MSIN, MCOS 

AUTHOR: JACQ 

SOURCE: SYM I 

OBJECT: Relocatable in Block "MATH" 



PURPOSE 

To find the sine or the cosine of a mid-precision floating point number and 
leave the result in the same format in the software registers MNT1, MNT2, 
MNT3. 

USAGE 



Calling Sequence 

The routine returns at L+2. 



L-1 


SMB 


MSIN (MCOS) 


L 


JSX 


MSIN (MCOS) 


L+l 


D 


ARG 


L+2 


Return 




Argument De 


scription 





The argument is given in radians, and is a mid-precision floating point 
number occupying three consecutive words in memory, beginning at address 
"ARG" as pointed out in the calling sequence. 

METHOD 



Let 

F (X) = sin X 

or F (X) = cos X 

Through transformation of the given argument X, the sine function of a pos- 
itive arc x is substituted to the given function F(X). The new arc x is de- 
fined in terms of two positive numbers, an integer N, and a fraction S. 
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cos X 


= sin(7T/2+ /Xj ) 


sin(-/x/) 


sinftT-X) 


M 


Q-7T/4 + S 


N 


Q -8K for sin X and X > 


N 


Q+4-8K for sin X and X <0 


N 


Q+2-8K for cos X 


X 


N-fT/4 + s 


F(X) 


= sin x 


0< N < 


8 


0< s < 


7774 


N defines a 


half-quadrant, from to 7. S is replaced by its complement 


if N is odd. 




Z 


= S if N is even 


z 


7772 -S if N is odd 


f(Z) 


sin Z if N = 0, 3, 4, 7 


f(Z) 


cos Z if N = 1,2, 5,6 


The function f(Z) is approximated by means of a 5-term Chebyshev poly- 


nomial: 




sin Z 
cos Z 


3 9 
S • Z + S • Z ... +S • Z + Se 
1 2 4 - 

C + C-Z . . .+C -Z + Ce 

o 1 4 — 


Finally the 


result is set negative, if N points to one of the last two quadrants: 


F(X) 


f (Z) if N < 4 


F(X) 


- -f (Z) if N > 4 


The coefficients of the polynomials are: 


Sine 


Cosine 


s o - ] 


C Q = .99999 99995 


S l = " 


.16666 6664 C = -.49999 9973 


S 2 = 


.83333 1376 E-2 C = .41666 4514 E-l 


S 3 = " 


.19835 7897 E-3 C = -.13882 8607 E-2 


S 4 = 


.26931 0165 E-5 C = .24112 6543 E-4 


sin Z 


" 31 Ce<2- 31 
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ERROR CONDITIONS 



13 
For argument values greater than or equal to 2 the function is set 

to zero, not computed. 

To form the function for arguments of this size, extra instructions 
would be required that would unfavorably affect the timing of the 
routine for all argument values. Since the results for large argu- 
ments would be of such poor accuracy (see following section) the 
function is arbitrarily set to zero. 

ACCURACY 



29 bits for arguments whose magnitude is under 1. Deduct 1 bit for each 
positive unit in the argument exponent, down to 16 bits. Beyond that arbit- 
rary limit, the function is replaced by zero. 

RESTRICTIONS 

Loading 

This subroutine locally references storage words and constants located in 
the "MATH POOL" module, and must therefore be loaded in the same 2K 
block as the pool. 

Other Routines 

DAD, DIVS, DSM, D2C.FLD, FSB, M.P, M. ZE, MPYS. 

Space Used 

119 words. 

Timing 

Hardware Multiply &: Divide 

1990 
1910 
2350 
2110 





Standard 


Minimum (MS IN) 


2980 


Minimum (MCOS) 


2730 


Maximum 


3800 


Average 


3110 
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SINE AND COSINE - MP FLOATING POINT 




K = K + Q 
£ Z <TT/4 



F(X) = 



I 



EXIT 



1 



z = x 



J 
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ODD 



QUADRANT s K/2 



EVEN 



I 



SUBROUTINE "M. P" 

SIN Z " il_ S n * z2n+1 
n*0 



I 



Y = SIN Z 




ODD 



1 



K x KS + 8i 
0-6KS^8 



z.I-z 

4 



QUADRANT = 



K + 1 



1 



/SUBROUTINE "M. P"\ 
n«l 



I 



Y = COS Z = C + C 




EXIT 
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Drawing No . 
390007 



ID Code 



BSJ 
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7GO PROGRAMMING SYSTEMS 



ARC TANGENT - MP FLOATING POINT 



QUALITY SOFTWARE 

DATE: April 1968 

ID CODE: BSM 

DRAWING: 390010 (Rev B) 

LABEL: MATN 

AUTHOR: j A CQ 

SOURCE: SYM I 

OBJECT: Relocatable in Block "MATH" 



PURPOSE 

To obtain the arc tangent of an argument given in the mid-precision floating 
point format. The result is set in the same format in the software registers 
MNT1, MNT2, MNT3. 

USAGE 



Calling 


Sequence 




L-1 


SMB 


MATN 


L 


JSX 


MATN 


L+l 


D 


ARG 


L+2 


Return 




Argument Description 





The argument is a mid-precision floating point number occupying three con- 
secutive words in memory, beginning at address "ARG" as pointed out in the 
calling sequence. 

METHOD 

To improve the timing with negative tangents, all computations are perform- 
ed on the absolute value of the arguments, then the final arc is given the sign 
of the original argument. Let 



Y 


= 


arctan X 


X 


= 


1x1 


X 

Z 


= 


1*1 

X 


Z 


= 


x-tan Ax 


1+x • tan Ax 



if |x| < 1 
if |x| > 1 

if x < 1/8 
if x > 1/8 



BSM - Continued Page 2 

Where tan Ax is a multiple of 1/8, such that: 

x-1/8 < tan Ax< x 
Then 

< Z < 1/8 

arctan x = arctan Z + Ax 

To obtain arctan Z, the Taylor series is expanded in terms of Z, trun- 
cated to the first five terms because of the reduced range of Z. 

The value of Ax is found in radian in a seven-entry table. 
Finally: 

y = arctan x if |X| < 1 

y = 1772 -arctan x if |x| > 1 

Y = y if X > 

Y = -y if X < 

ERROR CONDITIONS 
None 

RANGE AND ACCURACY 

The arc is computed in the interval -Tl/2, +1T/2. 
The accuracy is 29 bits. 

RESTRICTIONS 

Loading 

This subroutine locally references storage words and constants located in 
the "MATH POOL" module, and must therefore be loaded in the same 2K 
block as the pool. 

Other Routines 



ACMY, DLD, DSM, DST, D2C, FAD, FDV, FLD, FNRM, M.P, MPYS 



BSM - Continued Page 3 

Space Used 
128 words 
Timing (in cycles) 

Standard Hardware Multiply &: Divide 

Minimum 2700 2000 

Maximum 4800 2900 

Average 4120 2690 

tangent under 1 : 

Average 3820 2470 
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A IS A MULTIPLE 
OF 1/8 



TAN" A IS FOUND 
IN A 7 ENTRY 
TABLE 



Z 9 


T . A 
1 + A. 


X 


8 — 


^l 




o£ 


8 





I 



SUBROUTINE "M. P" " 

TAN" 1 Z ^(-1)" Z 2n+1 
n=0 2n+l 



TAN'j T = 

TAN Z + TAN A 



FLAG COTANGENT 
SUBSTITUTION 
T = 1/ | X | 







,-1 ^ 




NO . 


/'COTANGENT^ 
^v^SU BSTIT U TIONX 

to4 


YES 




Y = 




. -i 












IT/2- TAN "T 
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DATE: April 1968 

ID CODE: BSL 

DRAWING: 390009 (Rev D) 

LABEL: MEXP 

AUTHOR: JACQ 

SOURCE: SYM I 

OBJECT: Relocatable in Block "MATH" 



PURPOSE 

To find the exponential value of a mid-precision floating point number and 
leave the result in the same format in the software registers MNT1, MNT2, 
MNT3. 

USAGE 



Calling Sequence 

The routine returns at L+2. 



L-1 


SMB 


MEXP 


L 


JSX 


MEXP 


L+l 


D 


ARG 


L+2 


Return 




Argument Dt 


ascription 





The argument is a mid-precision floating point number occupying three con- 
secutive words in memory, beginning at address "ARG" as pointed out in the 
calling sequence. 

METHOD 

Let the magnitude x of the argument X be redefined in terms of Ln 2 by 
means of a positive integer N and a positive fraction Z. 

x - |x| = N • Ln 2 + Z 

< N 

< Z < Ln 2 

x N.Ln 2 Z N Z 
e = e . e = 2 . e 
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Let Z be split in terms of Az, 


a multiple of 1/8, and whatever remains in Z. 


Z = Az + z 




Az £ Z < Az + 1/8 





< z < 1/8 

Z Az z 

e = e • e 

Az z 

The value of e is found in a five -entry table. The value of e is approx- 
imated by means of the Taylor series truncated to the first six terms. 

ERROR CONDITIONS 

The subroutine M.OV is called when the argument is greater than 
+88.0296918. 

RANGE AND ACCURACY 

The argument should be within the range -89.4159863, +88.0296918. 

The accuracy is 29 bits. 

RESTRICTIONS 

Loading 

This subroutine locally references storage words and cons tants located in 
the "MATH POOL" module, and must therefore be loaded in the same 2K 
block as the pool. 

Other Routines 

DAD, DIVS, DSM, D2C.FAD, FDV, FLD, FMP, FNRM, FST, M.OV, M.P, M. ZE, 

MPYS. 



Space Used 
142 words 
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Timing (in cycles) 

Standard Hardware Multiply & Divide 

Minimum 50 50 

Maximum 4750 2760 

Average 3950 2530 
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( MEXP y 



I 



U ^ Q + R 
Ln2" 



P A MULTIPLE 
OF 1/8 



e IS FOUND 
IN A 6 ENTRY 
TABLE . 



FETCH ARGUMENT X 




u» |x| 




I 



R = P + Z 
5 

8 

8 



0& Z^. ■ 



SUBROUTINE "M.P" 



"TBI n! 



e U - (SMI e P - 2 Q 




.U 



EXIT 



1 



Q a 
R = U 



I 





OVERFLOW 

"U.OVi' 
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QUALITY SOFTWARE 

DATE: April 1968 

ID CODE: BSK 

DRAWING: 390008 

LABEL: MLOG 

AUTHOR: JACQ 

SOURCE: SYM I 

OBJECT: Relocatable in Block "MATH" 



PURPOSE 

To find the natural logarithm of a mid-precision floating point number and 
leave the result in the same format in the software registers MNT1, MNT2, 
MNT3. 

USAGE 



Calling Sequence 



The routine returns at L+2, 



L-l 


SMB 


MLOG 


L 


JSX 


MLOG 


L+l 


D 


ARG 


L+2 


Return 




Argument D 


escription 





The argument is a mid-precision floating point number occupying three con- 
secutive words in memory, beginning at address "ARG" as pointed out in the 
calling sequence. 

METHOD 

Let the argument be 

X = Xm ' 2 Xe 

.5< Xm < 1 



Ln X = Ln Xm + Xe • Ln 2 



BSK - Continued 


Page 2 


Let Xm be redefined as the product of two factors, Ax and x, 


Ax a multi- 


pie of 1/8: 




Xm = Ax • x 




1/2 < Ax < 7/8 




0< Xm-Ax < 1/8 





Ln Xm = Ln Ax + Ln X 



Ln Ax is obtained from a 4- entry table. Let 



z = 



Xm - Ax 
Xm + Ax 



0< z < 1/8 



x = 



1 + z 



Subtracting the two Taylor series of Ln(l+z) and Ln(l-z): 
Ln x = 2(z+z 3 /3 . . . +z 9 /9) + E 



E < 2 " 31 

ERROR CONDITIONS 

A negative or a null argument will result in a call to the subroutine M. OV. 

RANGE AND ACCURACY 

Any positive argument is acceptable. The result is accurate to 29 bits. 

RESTRICTIONS 

Loading 

This subroutine locally references storage words and constants located in 
the "MATH POOL" module, and must therefore be loaded in same 2K block 
as the pool. 



BSK - Continued p age 3 

Other Routines 

DDV, DSM, FAD, FLD, FLT, FMP, FSB, FST, M. OY, M.P 

Space Used 

83 words 

Time (in cycles ) 

Standard Hardware Multiply &c Divide 

Minimum 4000 2700 

Maximum 4600 2850 

Average 4230 2740 
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NATURAL LOGARITHM - MP FLOATING POINT 



I MLOG V- 



F A MULTIPLE 
OF 1/8 

.. * + Z 



Ln F IS FOUND 
IN A 4 ENTRY 
TABLE. 



FETCH ARGUMENT 
X • X • 2 X « 




X « F • U 
m 

2 8 



z = 



X +F 



o*z«I 



REVERSE SIGN 
OF ARGUMENT: 



/SUBROUTINE "M. P" 
\ n«0 2n+l 



SUBROUTINE "M.OV" 
TREAT EITHER CASE 
AS OVERFLOW 



J 


[ 




Ln X = 2 • S + Ln F 
m 

Ln X = Ln X + X . Ln2 
m e 

















EXIT 
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DATE: May 1968 

ID CODE: BSS 

DRAWING: 3 90016 

LABEL: M. P 

AUTHOR: J ACQ 

SOURCE: SYM I 

OBJECT: Relocatable in Block "MATH" 



PURPOSE 

To assist the mathematical function subroutines in the expansion of series 
in the mid-precision floating point format. 

This subroutine can sum up: 

2 . +An. X n 



P = X + A2-X< . .. .-_. 

or P* = Al.X 2 + A2.X 4 ...+An.X 2n 
or P = X + A1.X3 ...+An.X 2n+1 



USAGE 



The use of M.P is not encouraged outside the library. The reader is 
referred to the assembly listings of MEXP, MLOG, MSIN, MATN for 
examples of calling sequences and coefficient tables. 

Error Conditions 



Overflow conditions are flagged by FAD and FMP, routines called by M.P. 
RESTRICTIONS 

Loading 



M.P must be loaded in the same 2K block as the "MATH POOL". 



Other Routines 



FAD, FMP, FNRM, FST 

Space Used 
40 words 



BSS - Continued Page 2 

Timing (in cycles) 

Standard Hardware Hardware Multiply 

PI 398 + 507 (n-1) 233 + 342 (n-1) 

P2 307 + 507n 142 + 342n 

P3 291 + 507n 126 + 342n 

Deduct or add 20 cycles per term positive or negative. 
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RAYTHEON 



700 PROGRAMMING SYSTEMS 



DOUBLE SHIFT MAGNITUDE 



QUALITY SOFTWARE 



DATE: May 1968 

ID CODE: EST 

DRAWING: 390017 

LABEL: DSM 

AUTHOR: JACQ 

SOURCE: SYM I 

OBJECT: Relocatable in Block "MATH" 



PURPOSE 

To assist the mathematical function subroutines in the double arithmetic 
shift operations. 

USAGE 



The use of DSM is not encouraged outside the library. 

The shift is performed on the contents of the software registers MNT2, 
MNT3. The sign bit of MNT3 is left out of the shift. On return both 
sign bits are cleaned off and MNT2 is duplicated in the accumulator 
register. 

Calling Sequence 

Place the right byte of the shift instruction (left or right, double, shift 
count up to 15) in the accumulator register, before calling the routine: 

L-l SMB DSM 

L JSX DSM 

L+l Return 

Error Conditions 



No check is made on the right half of the shift instruction provided by the 
caller. 
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RESTRICTIONS 



Loading 

This subroutine locally references storage words located in the 
"MATH POOL" module, and must therefore be loaded in the same 2K 
block as the pool. 

Other Routines 



None 

Space Used 

12 words 

Timing ( in cycles) 



Maximum 26 
Minimum 23 
Average 24 
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700 PROGRAMMING SYSTEMS 



MP HYPERBOLIC TANGENT 



Quality software 



DATE: 

ID CODE: 

DRAWING: 

LABEL: 

AUTHOR: 

SOURCE: 

OBJECT: 



April 1969 
BYB 

392279 

MTNH 

JACQ 

SYM I 

Relocatable in block "Math" 



PURPOSE 

To find the hyperbolic tangent of a mid -precision floating point number and 
leave the result in the same format in the software registers MNT1, MNT2 
MNT3. 

USAGE 



Calling Se 


squence 




L-1 


SMB 


MTNH 


L 


JSX 


MTNH 


L+l 


D 


ARG 


L+2 


Return 




Argument 


Description 



The argument is a mid-precision floating point variable occupying three con- 
tiguous words in memory, beginning at address "ARG" as indicated in the 
calling sequence. 

METHOD 

Operations are carried on the absolute value X of the variable x, then the 
function is given the sign of x. 

The range of X = Ixl is divided into four intervals: 
, -,-29 



X 



tanh x = x 



•29 



^ X ^ 1 



tanh x = 



2x 
e -1 

e +1 
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2x 
1 ^ X ^.16 tanh x = € 



2x +i 
e +1 



16 ^E. X tanh x = — 

.X. 

In order to keep the accuracy of the function when X is within the interval 

-29 2X 

(2 , 1), it is necessary to obtain the difference e -1 directly, rather 

than subtracting 1 from e 

The exponential routine "MEXP" works as follows: 

2X = Q-Ln 2 + JF + Z (mod 1/8)] (mod Ln 2) 

2X F Z Q 

e = e • e • 2 

A special exit is provided in "MEXP" to give the partial result Y: 

Y Me Z -De F 

2x 
Upon return from MEXP the expression e -1 is obtained as follows: 



2X 
e 



> = [Y + ,e F -»] 2 Q + U Q -„ 

ERROR CONDITIONS 

None 

ACCURACY 

28 magnitude bits 

RESTRICTIONS 

Loading 

This subroutine locally references storage words and constants located in 
the "Math Pool" module, and must therefore be loaded in the same 2K block 
as the pool. 
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MEXP References 

Three non-entry words of the subroutine MEXP are referenced in MTNH 
as relocatable addresses. If these words are displaced in MEXP, their 
addresses should be redefined in MTNH. 

Other Routines 

D2C, FAD, FDV, FLD, FSB, FST, MEXP 

Space Used 

91 words 

Timing (Machine Cycles) 

Without Hardware With Hardware 

Multiply and Divide Multiply and Divide 

Minimum (X^2~ 9 ) 35 35 

Maximum 4800 3100 

Average 4550 2950 
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